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Pt nanocontacts, like those formed in mechanically controlled break junctions, are shown to de- 
velop spontaneous local magnetic order. Our density functional calculations predict that a robust 
local magnetic order exists in the atoms presenting low coordination, i. e., those forming the 
atom-sized neck. In contrast to previous work, we thus find that the electronic transport can be 
spin-polarized, although the net value of the conductance still agrees with available experimen- 
tal information. Experimental implications of the formation of this new type of nanomagnet are 
discussed. 
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Fabrication of metallic nanocontacts permits to probe 
the electronic and mechanical properties of conventional 
metals with unconventional atomic coordination^. Elec- 
tron transport in these systems depends on the tiny frac- 
tion of atoms in the sample forming the atom-sized neck 
which have a reduced coordination and are responsi- 
ble for the two-terminal resistance. Transport experi- 
ments can thereby probe the atomic and related elec- 
tronic structure of these atoms and provide informa- 
tion about a fundamental question: How bulk properties 
evolve when the system reaches atomic sizes and atoms 
with full bulk coordination are no longer majority. A 
bulk property that is susceptible to change is magnetism. 
Bulk Pt, for instance, is a paramagnetic Fermi liquid with 
a rather large spin susceptibilitj*^. A transition to a fer- 
romagnetic state could be expected upon reduction of 
the atomic coordination with the concomitant increase 
of the density of states at the Fermi energy beyond the 
Stoner limit. Density functional calculations^''*'^ for one- 
dimensional infinite Pt chains support this hypothetical 
scenario, resulting in a ferromagnetic transition above a 
critical lattice spacing which, depending on the compu- 
tational approach, can be below the equilibrium lattice 
constanljii. The formation of local moments in real Pt 
nanocontacts would not be totally unexpected. 

Formation of and electronic transport in finite Pt 
chains, for instance, have been extensively studied 
experimentalljiSiSiifliiiii^. Based upon the appearance 
of a peak at G = 0.5 x 2e^/h = O.SGq in the conduc- 
tance histogram Rodrigues et al. suggested that Pt and 
Pd nanocontacts could be spin polarizedpi^. The origin 
of this peak has been later attributed to adsorbateslSi 
so that magnetism in Pt and Pd nanocontacts has not 
been confirmed experimentally yet. Previous theory 
work has addressed the formation of local moments in 
Pd nanocontacts- and in Co , Pd and Rh short chains 
sandwiched between Cu planes^. To the best of our 
knowledge theory work on Pt nanocontactsi^*i^iiS4iLi& 
has overlooked the possibility of local magnetic order so 
far. Here we perform density functional calculations of 
both the electronic structure and transport and find that 
local magnetic order can develop spontaneously in Pt 



nanocontacts. Local magnetic moments as high as 1.2 fiB 
in low-coordination atoms are found. Interestingly, while 
transport is definitely spin polarized, the calculated total 
conductance of magnetic and non-magnetic Pt nanocon- 
tacts is very similar and in agreement with experimental 
data, explaining why magnetism has been unnoticed so 
far. 

The rest of the paper is organized as follows. First we 
discuss the methodology of the electronic structure cal- 
culation presented here. Second we present results for 
idealized infinite and finite one-dimensional Pt chains. 
Then we calculate both zero-bias transport and the elec- 
tronic structure of Pt nanocontacts, which turn out to be 
magnetic beyond attainable values of the stretching. Fi- 
nally we discuss whether magnetism in Pt nanocontacts 
is robust to thermal and quantum fluctuations, and pro- 
pose experimental verification of the different scenarios. 
We also comment on the effect of spin-orbit interaction, 
missing in our calculations. 

Electronic structure and transport calculations.- The 
electronic structure of various low dimensional struc- 
tures of Pt which mimic actual nanocontacts are cal- 
culated in the density functional approximation, using 
either CRYSTALOSia or ALACANT (ALicante Ab ini- 
tio Computation Applied to Nanotransport), our ah ini- 
tio transport package^" interfaced to GAUSSIANOSSi. 
ALACANT describes the bulk electrodes using a semiem- 
pirical tight-binding Hamiltonian on a Bethe lattice. 
The spin-resolved density matrix includes the effect of 
the electrodes by means of self-energies and the spin- 
polarized transportSSiS^ is obtained with the standard 
Landauer formula. Both GAUSSIAN03/ ALACANT and 
CRYSTAL03 perform electronic structure calculations 
using a basis of localized atomic orbitals (LAO). CRYS- 
TAL03 permits to calculate infinite systems with crys- 
talline symmetry whereas ALACANT is suitable for sys- 
tems with no symmetry or periodicity such as nanocon- 
tacts in or out-of-equilibrium situations. In the case of 
the Pt calculations shown here we use scalar relativis- 
tic (SR) pseudopotentials for the 60 inner electrons and 
the remaining 18 electrons are treated using generalized 
gradient approximation (GGA) density functionals. The 
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FIG. 1; (Color online), (a) Energy per atom for a perfect 
monostrand Pt chain, (b) Magnetic moment per atom as a 
function of lattice spacing a. (c) Density of states (D.O.S.) of 
the paramagnetic chain at the Fermi energy as a function of 
a. (d) D.O.S. as a function of energy for a — 2.4A (dashed) 
and a = 2.8A (solid) 



basis set used for all the calculations has been optimized 
to describe bulk Pt as well as Pt surfaces^. Other ba- 
sis sets such as LANL2DZ or SDD— have occasionally 
been employed for comparison. The main results do not 
depend on the choice of basis set. 

Results.- We first consider a perfect one-dimensional 
mono-strand Pt chain. Such an idealized system serves 
as a standard starting point to understand lower symme- 
try geometries. It also permits to test whether our LAO 
pseudopotential methodology reproduces the results ob- 
tained with SR all electron plane-wave calculations re- 
ported by Delin et alA. In Fig. ^a) we show the energy 
per atom as a function of the lattice constant a both 
for the paramagnetic (PM) and the ferromagnetic (FM) 
chain. They both have a minimum at a = 2.4A. The FM 
chain develops a non-negligible magnetic moment when 
the lattice constant goes beyond a ~ 2.6A. This configu- 
ration is clearly lower in energy above that distance. The 
energy difference between the FM and the PM configu- 
rations is 16 meV per atom for a = 2.7A and 33 meV 
per atom for a = 2.8A. The magnetic moment per atom 
reaches a saturation value of 1.2/xb. The equilibrium dis- 
tance, critical spacing, asymptotic magnetic moment and 
shape of the phase boundary obtained by us are similar 
to those obtained by Delin et al.^ using a SR all-electron 
plane-wave calculation. Our results and those of Delin 
et al. underestimate the onset of the magnetic transition 
compared to calculations including spin-orbit coupling^'^ 
that predict that a magnetic moment forms already be- 
low the equilibrium distance. 

The magnetic transition in the phase diagram [Fig. 
1(b)] is compatible with the Stoner criterion for ferro- 
magnetic instability. As the chain is stretched, the atom- 
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FIG. 2: (Color online). Energy bands for ideal Pt chain with 
a = 2.8A. Left: ferromagnetic phase. Right: paramagnetic 
case. 



atom coupling becomes weaker, the bands narrow down 
and so does the density of states (p(e)). Since the in- 
tegrated p(e) must be equal to the number of electrons 
per atom, narrowing of the DOS implies an increase of 
the p(e) [see Fig. ^d)] and, therefore, an increase of the 
spin susceptibility, which is proportional to pitp)- In Fig- 
n^c) we show how the pitp) of the PM chain increases 
as a function of the lattice constant. The remarkable 
feature of Pt chains is that the Stoner instability occurs 
close to the equilibrium lattice spacing. 

The electronic structure of the ideal Pt chain sheds 
some light on the electronic structure of the nanocon- 
tact. In Fig. |21we show the energy bands for the ideal Pt 
chain both in the FM (left panel) and PM (right panel) 
configurations, for a lattice spacing a = 2.8A. We notice 
that the spectrum at /c = has four resolved energy lev- 
els per spin. These correspond to the 6s level, and the hd 
levels which, because of the axial potential created by the 
neighboring atoms, split into 2 doublets Ei and one 
singlet A\. The E\ and Ei are linear combinations of or- 
bitals with = ±1 and = ±2, respectively, whereas 
the A\ singlet is a = orbital that hybridizes with 
the lower energy 6s orbital. The largest contribution to 
the density of states, and therefore to the magnetic in- 
stability, comes from the >li-like band at the edge of the 
Brillouin zone. However, the prominent role played by 
these bands in the magnetic behavior of Pt chains is in 
stark contrast with their role on the transport proper- 
ties of Pt nanocontacts (see below and see also related 
work on Ni nanocontacts^^) . Four spin-degenerate bands 
cross the Fermi energy in the PM case whereas 7 spin- 
split bands do it in the FM chain. In the FM chains the 
number of spin minority channels is 6, and the number 
of spin majority channels is 1. Although spin-orbit in- 
teraction modifies significantly the bands*, the number 
of bands at the Fermi energy is pretty similar in both 
cases. Therefore, one can anticipate that the number 
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of open channels in the magnetic and non-magnetic Pt 
nanocontacts studied below should be roughly the same 
and thereby the conductance should be similar, although 
the spin polarization might well be large in the former 
case. The conductance of the ideal FM chain is 3.5Gq, 
very far from the value of O.SGq that allegedly signals the 
emergence of magnetism and also far away from half of 
the conductance of the PM chain, so it is very unlikely 
that the celebrated half quantum can be attributed ex- 
clusively to magnetism. 

Real Pt chains are typically less than five atoms long 
and are connected to bulk electrodes. Although not sur- 
prising, we have verified that magnetism survives in iso- 
lated short chains with iV^ = 3, 4 and 5 Pt atoms. The 
equilibrium distance is 2.4A for all Na = 3, — 4, 
and Na — 5. Interestingly, the short chains are always 
magnetic in the Na — 3 and Na — 4 cases and show a 
non-magnetic to magnetic crossover at a = 2.6A in the 
Na = 5 case, already similar to the ideal infinite chain. 
The total magnetic moment of all the Na = 3 chains 
with a < 3.0A is 4/iB. The outer atoms have a magnetic 
moment of 1.36/iB and the central atom with larger co- 
ordination has a smaller magnetic moment of 1.29/iB. In 
the case of Na = 4 the total magnetic moment is 6/iB, 
and their distribution is similar to the Na = 3 case. 

The calculations above show that magnetism is present 
both in finite- and infinite-sized Pt systems with small 
atomic coordination. It remains to be seen that this 
holds true in nanocontacts where none or only few atoms 
have a small coordination (like in the case of formation of 
short chains) , but these are strongly coupled to the bulk. 
In order to verify whether or not this is the case, we 
have calculated both electronic structure and transport 
for a model Pt nanocontact. It is formed by two opposite 
pyramids grown in the (001) crystallographic orientation 
of bulk Pt and joined by one atom which presents the 
lowest possible coordination [see inset in Fig. OJa)]. Re- 
laxation of the 11 inner atoms of the cluster has been 
performed starting from an equilibrium situation as a 
function of the distance d of the outer planes. Zig-zag 
configurations appear in the chain for small values of d 
(not shown until the three-atom chain straightens up 
(see inset in Fig. OJa)] followed by a plastic deformation 
(not shown). This deformation can be in the form of a 
rupture or a precursor of the addition of a new atom to 
the chain which comes from one of the two 4-atom basest. 

We now compute the transmission before the plastic 
deformation occurs (left panel in Fig. |3J, where the 
atom-atom distance in the short chain is 2.82 A. The 
value of the corresponding atomic-plane-averaged mag- 
netic moments are also shown in the inset. As expected, 
it decreases for atoms in the bulk as the coordination 
reaches the bulk value. Some atomic realizations in the 
stretching process (like zig-zag ones) result in nanocon- 
tacts with smaller Pt-Pt distance and no magnetism, in 
agreement with the infinite chain phase diagram in Fig. 
1. In contrast to the infinite chain, there are only three 
channels contributing to the total conductance for mi- 
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FIG. 3: Conductance per spin channel for a nanocontact 
with a 3-atom Pt chain (see inset) for the magnetic solution 
(a) and the non-magnetic one (b). The atom-atom distance 
in the chain is 2.8 A. 



nority electrons and there are more than one (three) for 
majority ones. For the majority electrons these are a 
perfectly transmitting s-type channel and two partially 
open (T = 0.4) pd-type channels (one Pxdxz- and one 
Pydyz- hybridized). The three minority channels have the 
same character as the majority channels, except that here 
the s-type channel does not transmit perfectly while the 
transmission of two pd-type channels is enhanced so that 
all three minority channels have a transmission around 
0.7. The other two remaining pd-like channels are respon- 
sible for the sharp resonances that appear around the 
Fermi level. The total conductance of the nanocontact 
in the FM case thus turns out to be around 4e^//i = 2Go 
which is only slightly larger than the average experi- 
mental value corresponding to the last plateau (I.TSGq), 
but, interestingly, barely differs from the value obtained 
when the possibility of magnetic order is ignored [around 
2.3Go, see right panel in Fig. As in the case of the 
ideal chain, the FM conductance is not half of the PM 
conductance nor half of Go- To conclude this discussion 
we notice the, although transport is only weakly spin 
polarized, magnetism brings the pd-\ike resonances up 
to the Fermi level compared to the non-magnetic case. 
These resonances may well give features in the low bias 
conductance not present if Pt were not magnetic. 

Conclusions and Discussion.- The main conclusion of 
this work is that density functional calculations pre- 
dict that nanochains formed in Pt nanocontacts can be 
stretched as to become magnetic. The magnetic moment 
is localized mainly in the atoms with small coordination 
and does not modify appreciably the total conductance, 
although the transmission is moderately spin polarized. 
How robust are these results? It is well known that both 
local and gradient-corrected density functionals present 
some degree of electronic self-interaction, in contrast 
with the Hartree-Fock approximation. Self-interaction 
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is larger for localized electrons and shifts the the d bands 
upwards in energy, as shown in the case of Co, Ni and Pd 
one dimensional chains^^ . A number of schemes to avoid 
this problem, like LDA+U and Self-Interaction Correc- 
tion functionals have been proposed. The method of 
choice between chemist is hybrid functionals^^ in which 
local and Hartree-Fock exchanged are combined and the 
self-interaction is reduced. We have calculated the mag- 
netic phase diagram of the one dimensional Pt chain us- 
ing the hybrid B3LYP functionalS& and found, somewhat 
expectedly, that magnetism is enhanced and that B3LYP 
infinite Pt chains are ferromagnetic down to the equi- 
librium distance (a = 2.4 A). Both non-local exchange 
and spin-orbit coupling^ enhance the stability of mag- 
netism in Pt nanocontacts. This and previous results on 
Ni nanocontacts^'^ lead us to believe that self-interaction 
is an issue in the electronic structure and transport prop- 
erties of transition metal nanocontacts and further work 
is necessary along these lines^. 

The mean field picture of the electronic structure de- 
scribes a static magnetic moment without preferred spa- 
tial direction. In reality the nanomagnet formed in the 
break junction is exchanged coupled dynamically to the 
Fermi sea of the conduction electrons of the electrodes. 
In the case of a spin S = 1/2 this can result in the for- 
mation of a Kondo singlet that would yield an anomaly 
in the zero bias conductance. For larger spins, the con- 
duction electron sea cannot screen the spin completely 
so that the magnetic moment survives. The magnetic 
moment of the nanocontact in Fig. 3 is 5 w 6, compa- 



rable to that of single molecule magnetsSS, making the 
formation of a Kondo singlet unlikely. 

In the absence of spin-orbit interactions and external 
magnetic field a electronic configuration with total spin 
S has 25 -|~ 1 degenerate configurations corresponding to 
the spin pointing along different directions. However, 
spin-orbit interaction is strong in Pt and produces spin 
anisotropy, favoring orientation along the transport di- 
rection axis in the case of one dimensional chains^. Ther- 
mal fluctuations of the magnetic moment between these 
two configurations are quenched for temperatures smaller 
than the anisotropy barrier. Departures from the easy 
axis orientation will be damped via electron-hole pair 
creation across the Fermi energy that would also result 
in small bias features in transport^ ^. The application of 
a sufficiently strong magnetic field in the direction per- 
pendicular to the easy axis moves the local magnetic mo- 
ments away from their easy axis. This is known to change 
the number of open channels at the Fermi energy in both 
the case of Ni ideal chainsSS and in the case of ferromag- 
netic semiconductor tunnel junctions^. This effect, or 
maybe even larger, can be expected in Pt nanocontacts 
and could be used to detect the nanomagnetism experi- 
mentally. 
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